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Figure 1. BFACF moves on a part of a square or simple cubic lattice polygon. 
Moves of type I are positive elementary moves when 2 edges are added, and 
negative elementary moves when 2 edges are removed. Moves of type II do not 
change the length of the polygon, and are called neutral elementary moves. 



1. Introduction 

The BFACF algorithm [1, 2, 5] is a Metropolis style Monte Carlo algorithm [17] that 
samples self-avoiding walks with fixed endpoints in the grand canonical ensemble in 
the hypercubic lattice Z d . The algoritm uses the local elementary moves in Figure 1 
to sample walks along a Markov Chain from the (non-Boltzman) distribution 

P ne ' 5n (1) 

" ]>>c„(x, y ) e ^ 

n>0 

where (3 is a parameter and c n (x, y) is the total number of self-avoiding walks of length 
n from the lattice site x to the lattice site y. 

The BFACF algorithm has also been used to sample unrooted polygons in the 
square and cubic lattices (see for example references [3, 4, 11, 12, 20, 21] and it has 
been generalised to lattice ribbons [19]. 

The BFACF algorithm can be used to sample knotted lattice polygons of fixed 
knot type [16] (also see [10] for a review) and has recently been used to determine the 
entropy and length of minimal knots in the cubic lattice [24, 25]. 

Despite its widespread use, the BFACF algorithm appears to have only been used 
in simple hypercubic lattices . In this paper our aim is to extend the elementary moves 
of the BFACF algorithm (see Figure 1) in the simple cubic lattice (hereafter referred 
to as the SC) to the body-centred and face-centred cubic lattices (hereafter referred 
to as the BCC and FCC lattices respectively). In addition, we examine the ergodicity 
properties of the proposed elementary moves when applied to unrooted self-avoiding 
polygons in the BCC and FCC lattices. 

The BFACF algorithm is known to have elementary moves with non-trivial 
ergodicity properties in the cubic lattice [16, 8]. In this paper we prove an analogous 
result for the BCC and FCC lattices. 

BFACF moves on walks or polygons (see Figure 1) have been generalised by re- 
intcrprctation as plaquette atmospheres (see [13]) — namely the ways in which edges 
can be added, deleted or shuffled around plaquettef adjacent to edges in the polygon. 
Examining the set of such possible moves has proved extremely useful (see [13]) and 
inspired generalisations of the Rosenbluth algorithm [23] to the GARM [22] and GAS 
algorithms [14, 15]. 

f Unit squares in the lattice bounded by four lattice edges 
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There are natural analogues of these plaquette atmospheres on the BCC and FCC 
lattices and the new elementary moves follow immediately from the definitions of these 
atmospheres; see Sections 3.1 and 4.1. While these moves are quite easy to define, 
some work is required to demonstrate their ergodicity properties. 

In reference [15] an implementation of the BFACF elementary moves using the 
GAS algorithm was used to determine the shortest knots of given type in the simple 
cubic lattice. In this paper we extend those results by adding data for the BCC 
and FCC lattices, implementing the new BFACF-stylc elementary moves using the 
GAS algorithm. The minimal lengths of a selection of knot types in these lattices are 
displayed in Table 1. Knot types are indicated in the first column, while the next 
three columns give the lengths of the shortest knots of given type in each lattice. The 
last three columns are the number of distinct polygons (or "population" ) of the lattice 
knots of minimal length. 



Table 1. Minimal Knots in Cubic Lattices 
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In Section 2 we review the ergodicity properties of the BFACF elementary moves 
in the simple cubic lattice. We recall that the irrcducibility classes of these moves, 
when applied to unrooted simple cubic lattice polygons, are the knot types of the 
polygons as embeddings of the circle in R 3 . This result, in particular, implies that 
two unrooted cubic lattice polygons u>i and W2 are in the same knot type if and only 
if there is a sequence of (reversible) BFACF elementary moves which will take u>\ to 
0J2 [16]. We generalise this result to the BCC and the FCC lattices. 

In Section 3 we define an analogous set of elementary moves for polygons in the 
BCC lattice. We examine the properties of these moves and prove that the projection 
of any (unrooted) BCC polygon can subdivided by application of the BCC elementary 
moves. A corollary of this result is that polygons can be made contact free — roughly 
speaking the polygon can be "inflated" so that non-consecutive vertices do not lie close 
to each other. This result is then used to sweep the BCC polygon into a sublatticc 
isomorphic to the simple cubic lattice. In this sublatticc the BCC elementary moves 
reduce to the usual BFACF moves. This is sufficient to show that the BCC elementary 
moves are irreducible on the classes of unrooted BCC polygons of given knot type. 

Polygons in the FCC lattice are examined in Section 4. In this lattice the 
analogous of the cubic lattice BFACF moves is a single reversible elementary move 
which increases or decreases the length of an FCC lattice polygon by one. The proof 
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proceeds similarly to the BCC case, albeit with more cases. Again we show that 
any FCC lattice polygon can be made contact free and be swept onto a sublatticc 
isomorphic to the simple cubic lattice. The usual set of simple cubic lattice BFACF 
moves can be performed on the polygon (in this sublattice) by compositions of the 
FCC elementary move. Similarly, this is sufficient to show that the FCC elementary 
moves are irreducible on the classes of unrooted FCC polygons of given knot type. 

In Section 5 we show how the elementary moves may be used in a BFACF-style 
algorithm. This implementation samples from the distribution of polygons similar 
to that in equation (1). In addition, we note that the elementary moves can be 
implemented using the Metropolis algorithm instead, or using GAS-style sampling. 
We conclude the paper with a few final comments. 

2. Knotted Lattice Polygons 

Let S be the circle; we consider an injective map / : S — > R 3 , to be an embedding of 
the circle in Euclidean three space (that is, / is an injection, and is a homeomorphism 
onto its image). A polygon in the simple cubic lattice (SC), or the BCC lattice, or the 
FCC lattice, is a piecewise linear embedding of S into E 3 . Any such embedding of S 
is a knot, and if the embedding is a lattice polygon, then the embedding is a lattice 
knot. In this way cubic lattice polygons are lattice knots. 

Two oriented embeddings / and g are ambient isotopic if there is an orientation- 
preserving isotopy H : E 3 x I -> E 3 x I (where I = [0, 1]) with H(y,t) = (h t (y),t) 
such that ho is the identity (ho of — f) and the composition hi o / = g. In other 
words, two lattice polygons are ambient isotopic if there is a continuous deformation 
of E 3 which takes the embedding / of the first polygon onto the embedding g of the 
second polygon. 

Two lattice polygons in any cubic lattice (ie the SC, BCC or FCC lattices) are 
said to be equivalent if they are ambient isotopic. These equivalence classes of oriented 
embeddings of the circle into the cubic lattices define the knot types of lattice knots, 
see for example references [9] for a reviews and definitions of lattice knots. 

In this paper we prove that there exists piecewise linear realisations of orientation- 
preserving ambient isotopies between lattice knots of the same knot type in the 
BCC and FCC lattices. These isotopies can be constructed as sequences of local 
deformations of the lattice knots in terms of BFACF-style elementary moves. This is 
an extension of a similar result for polygons on the simple cubic lattice in [16]. In 
that result, the piecewise linear orientation-preserving isotopies are realised in steps 
using the elementary moves of the BFACF algorithm, illustrated in Figure 1. By 
constructing the isotopies in this way in reference [16] the following result is proven: 

Theorem 1 (Thm 3.11 from [16]) The irreducibility classes of the BFACF algo- 
rithm, when applied to unrooted polygons [in the simple cubic lattice], are the knot 
types of the polygons as piecewise linear embeddings in E 3 . <0 

This, in particular, follows by proving that for every pair of (oriented) polygons 
(ui, v) of the same knot type, there exists a finite sequence of elementary moves of 
either type I or II in Figure 1, such that these moves change oj into v. Each move is 
a local deformation of the polygon and the ambient space around it, and is itself an 
isotopy. The sequence of moves is a composition of these local isotopies, and is itself 
a realisation of an isotopy H : E 3 x / — > E 3 x I such that H(y,t) = (h t (y),t) where 
ho is the identity and h\L0 = v. 
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In other words, the equivalence classes of polygons induced by the simple cubic 
lattice BFACF elementary moves in Figure 1 coincides with the knot types of the 
polygons, proving that the irreducibility classes of the BFACF algorithm are the knot 
types of unrooted polygons in the simple cubic lattice. In this paper, we extend this 
result to the BCC and FCC lattices — these extensions are Theorems 6 and 10 below. 

2.1. Lattice Knot Projections 

The BCC and FCC lattices are generated by finite sets of basis vectors {cj} and {ej}. 
Once the origin O of the lattice is set, then vertices are defined by linear combinations 
^2iPiCi and J2iPi e i respectively in the BCC and FCC, where the Pi £ Z are finite 
signed integers. Thus, in each case the vertices have integer Cartesian coordinates. 

Two vertices u and v are adjacent if the difference u — v is a basis vector of the 
lattice. In this case an edge uv is defined between the vertices. The vertices u and v 
are the end-vertices of the edge uv. We consider a lattice to be the collection of all 
its vertices and edges. Two vertices are adjacent if they are the endpoints of the same 
edge. Two edges are adjacent (or incident on one another) if they share exactly one 
end- vertex. 

A lattice polygon is defined as a sequence or list of n adjacent edges 
(uoUi, U1U2, . . . , u„_iUo), such that all vertices {u , u 1; . . . , u„_!} are distinct. 

The length of the polygon is the number of edges n it contains (but its geometric 
length will generally be different from this, since the edges do not necessarily have 
length equal to one). 

The BCC lattice is Eulerian with girth 4, and all lattice polygons in it have even 
length. The FCC lattice is Eulerian of girth 3, and polygons of length longer or equal 
to 3 can be realised in this lattice. 

A lattice edge UjUj in the BCC lattice is said to be parallel to the Cj direction, 
or in the Cj direction if Uj — Uj = Cj. Similarly, one may define edges to be in the 
ej direction in the FCC lattice. When a lattice edge in the BCC is parallel to the 
Cj direction (or to the ej direction in the FCC), then we shall frequently abuse our 
notation by denoting it by its direction Cj (or by ej). 

A line segment in a lattice polygon is a maximal non-empty sequence of adjacent 
edges in the polygon of the form CjCj . . . Cj (in the BCC lattice), or e^e^ . . . ej (in the 
FCC lattice). We say that these line segments are in the Cj of e, directions respectively. 

In what follows, we shall work with the projections of lattice polygons winto 
geometric planes A along a direction u. To define these projections, consider two 
independent (unit) vectors co-planar with A, and let z = x x y be a vector normal to 
A. A vector u is transverse to A if u • z 7^ 0. 

The three vectors {x, y, u} is the basis of a (non-orthogonal) coordinate system S 
in R 3 . Points in the polygon uj can be identified by their coordinates in S, for example, 
uj is a piecewise linear curve parametrised by t and each point uj(t) has coordinates 
(x t ,y t ,u t ). 

The projection of uj into A along u is defined by the set of points (x t , y t ) in A for 
all values of the parameter t. 

A multiple point in the projection of uj into A along u is a point in the projection 
which is the image two or more distinct points in uj. A multiple point is a double point 
if it is the image of exactly two points. 

In the case of lattice polygons, projections of polygons will be subgraphs of the 
projection of the lattice into a plane normal to a given lattice axes. For example, the 
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Figure 2. Nine of the 12 plaquettes adjacent to an edge in the BCC lattice. Note 
that the six on the right are planar, while the three on the left are non-planar. The 
remaining three (undisplayed) polygons are mirror-images of the three non-planar 
polygons. 

projection of a simple cubic lattice polygon along the Z-direction into the XY-plane 
is a square lattice, and a cubic lattice polygon will project into a subgraph of this 
square lattice. This projection is a lattice knot projection in the square lattice; see 
reference [16]. 

In the BCC and FCC lattices we shall take projections of lattice polygons 
onto symmetry planes of the lattice, along directions which are transverse but not 
necessarily orthogonal to the symmetry planes. 

3. BFACF-Style Elementary Moves in the BCC Lattice 

In this section, we propose the local elementary moves of a BFACF-style algorithm in 
the BCC lattice and we show that they are sufficient to realise a piecewise linear 
orientation preserving isotopy on unrooted BCC polygons of the same knot type 
embedded in R 3 . In particular, this implies that the irreducibility classes of the 
BCC elementary moves coincide with the knot types of unrooted BCC polygons as 
determined by their embeddings in three space. 

3.1. BFACF style moves in the BCC lattice 

In section 2.1 the notion of lattice polygons and projections of polygons were defined 
in general. We note in particular that the basis vectors of the BCC lattice are points 
in ]R 3 with Cartesian coordinates given by pci + qc 2 + rc 3 + 5C4 where p, q,r,s G Z, 
and where the vectors Cj are given by 

Ci = (1,1,1), C 2 = (1,1,-1), C 3 = (1,-1,1), C 4 = (1,-1,-1), 

c 5 = (-1,-1,-1), c 6 = (-1,-1,1), c 7 = (-1,1,-1), c 8 = (-1,1,1). 

Observe that C5 = — Ci, eg = — C2, C7 = — C3 and eg = — C4. The vectors Cj the 
generating or basis vectors of the BCC and they all have (geometric) length 
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Figure 3. A projection of the BCC lattice and a BCC lattice polygon into the 
XY-plane. The projection of the BCC lattice is a square lattice rotated at 45° 
degrees with the Cartesian axes, and with edge-length \f2. The polygon projects 
to a subgraph of the projected lattice. 



Two vertices, u and v, are adjacent if they are separated by a single basis vector 
Cj. Lattice edges will be in the Ci directions in the BCC, and each lattice edge Cj will 
be an edge in minimal length polygons of length 4 in the BCC lattice. These minimal 
length polygons containing Cj bound 12 faces in the BCC, and these faces are not 
square and may not be planar in 6 of the cases. Nevertheless, we shall refer to them 
as lattice plaquettes - see Figure 2. 

The (orthogonal) projection of the BCC lattice into the XF-plane is a geometric 
square lattice with edge lengths \/2 and rotated at 45° with respect to the X-axis. A 
polygon in the BCC projects to a subgraph of this square lattice. This is illustrated 
in Figure 3. 

We define elementary moves on polygons in the BCC lattice by considering the 
12 plaquettes adjacent to every edge Cj (see Figure 2). Collectively, these plaquettes 
composed the plaquette atmosphere of the polygon [13]. In the simple cubic lattice, 
each edge is incident on at most 4 atmospheric plaquettes, and elementary moves of 
the BFACF algorithm are obtained by selecting an atmospheric plaquette and then 
exchanging edges along it boundary to update the polygon. 

In particular, by taking the alternate path around the boundary of atmospheric 
plaquette BFACF moves in Figure 1 is obtained. A type I move on the simple cubic 
lattice performed by replacing a single edge incident on an atmospheric plaquette P 
in the polygon by the 3 edges in the alternative path around the boundary of P. 
This move is reversible, and together the move and its reverse defines type I moves. 
Similarly, a type II move on the cubic lattice replaces a pair of edges on incident on 
an atmospheric plaquette P with the other two edges in P. 

The 12 atmospheric plaquettes incident on edges in BCC polygons will similarly 
be used to determine the set of elementary moves on BCC lattice polygons. A similar 
analysis to the above gives the following set of elementary moves on polygons in the 
BCC lattice: 

BCC Elementary Moves: 

• Replace an edge ci by c 2 cic 6 (and all permutations of these vectors in the set 
of vectors Cj). This move increases the length of the polygon by 2. Conversely, 
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Figure 4. Replacing ci by three edges as shown defines an elementary move 
which increases the length of a polygon in the BCC lattice by two. This defines 
a positive plaquette atmospheric move on the polygon. Reversing the move gives 
a negative plaquette atmospheric move which reduces the length of the polygon 
by two edges by replacing three edges on the right by the single edge on the 
left. All the possible permutations of the vectors gives the complete collection 
of positive and negative plaquette atmospheric moves. Observe that closing the 
vectors on the right with the dotted line gives a square in case I, and a non-planar 
quadrilateral in case II. 

replace C2C1C6 by Ci (and all permutations of these vectors) to obtain a move 
reducing the length of the polygon by 2. This moves is illustrated (generically) 
by case I in Figure 4. Observe that the vectors {c 2 ,Ci,c 6 } are coplanar, so 
that Cic 2 cic 6 forms a planar quadrilateral — a square. Hence, we call this the 
planar BCC positive and negative plaquette atmospheric moves, or the planar 
BCC positive and negative elementary moves. 

• Replace an edge ci by c 2 c 3 c 8 (and all permutations of these vectors). This move 
increases the length of the polygon by 2. Conversely, replace c 2 c 3 c 8 by Ci (and 
all permutations of these vectors) to obtain a move reducing the length of the 
polygon by 2. This move is illustrated by case II in Figure 4. Observe that the 
set of vectors {c 2 , c 3 , c 8 } is not coplanar; hence closing them off with Ci bounds 
a non-planar quadrilateral or atmospheric plaquette. These are the non-planar 
BCC positive and negative plaquette atmospheric moves, or the non-planar BCC 
positive and negative elementary moves. 

• Replace CiC6 by c@Ci (and all other permutations of these vectors). This is a 
neutral move, illustrated as case I in Figure 5. Note that these two pairs of 
edges form a square and so this move occurs in a plane in three space. These are 
the planar BCC neutral plaquette atmospheric moves, or the planar BCC neutral 
elementary moves. 

• Replace Cic 6 by c 3 c 2 (all other permutations of these vectors). This is a neutral 
move, illustrated as case II in Figure 5. Now, these two pairs of edges form a non- 
planar quadrilateral or atmospheric plaquette and so this move does not occur 
in a plane in three space. More generally, these are the non-planar BCC neutral 
plaquette atmospheric moves, or the non-planar BCC neutral elementary moves. 

Projecting the moves of Figures 4 and 5 into the A"F-plane gives the (two 
dimensional) square lattice BFACF moves shown in Figure 1, but in the rotated square 
lattice of Figure 3 instead. Thus, by executing the BCC positive, negative and neutral 
plaquette atmospheric moves of Figures 4 and 5 on a polygon in the BCC, the image 
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Figure 5. Replacing cgci by either I: cicg or by II: C3C8 defines neutral plaqucttc 
atmospheric moves. Observe that ci +C6 = C3 + cg. All the other possible neutral 
moves are obtained by using all the possible permutations of the in the above. 
Observe that by closing the vectors on the right into a quadrilateral along the 
dotted lines (which follows the path of cgci between the bullets) gives a square 
in case I and a non-planar quadrilateral in case II. 

of these moves in the projected square lattice include, as a subset, all the (usual simple 
square lattice) BFACF moves on the projected polygon. 

3.2. Stretching polygons in the BCC lattice 

We first outline our approach before we prove our main results. The proof of our main 
BCC lattice result in Theorem 6, is presented in two parts. We first show that any 
BCC lattice polygon can be swept onto a sublattice of the BCC which is isotopic to 
the simple cubic lattice (as an oriented embedded graph in three space). Then we 
show that in this sublattice we use the ergodicity properties of the simple cubic lattice 
BFACF algorithm [16] to complete the proof (see Theorem 1). 

In other words, we shall show that any BCC lattice polygon can be swept into 
a sublattice L with basis vectors {ci, c 3 , C4, c 5 , C7, c 8 }, and then show that a subset 
of the BCC elementary moves simulates a simple cubic lattice BFACF algorithm in 
this sublattice (which is not orthogonal but is nevertheless isotopic to the simple cubic 
lattice). 

Our approach would be to demonstrate that the BCC elementary moves are 
sufficient to replace polygon edges outside the sublattice L with edges in L, while 
avoiding any self-intersections in the polygon as it is updated in this process. This 
is achieved by stretching the polygon to create sufficient space for executing BCC 
elementary moves (sec Figure G). 

The stretching of a polygon will proceed by identifying a maximal line which 
intersect the projected polygon in its right-most and top-most projected edges. The 
polygon will be recursively stretched in directions normal to the maximal line by 
stretching parts of it across the maximal, while inserting edges to maintain its 
connectivity. We show that this can be done using the BCC elementary moves. 

In particular, project the BCC lattice along the Z-direction onto the XF-plane, 
and let w be a BCC polygon with projection Puj in the XF-planc. Then the image of 
the lattice and the polygon is the square lattice and a graph embedded in the square 
lattice illustrated in Figure 3. The maximal line K of Puj is the line x + y — k with k 
the maximum value such that K has a non-empty intersection with Puj. 
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Figure 6. Translating a line segment with a single positive move and then a 
sequence of neutral moves in the BCC lattice, projected into the Xy-plane. Since 
the projection of the polygon is to a square lattice, the BCC elementary moves 
project to square lattice BFACF elementary moves. 

K is the image of a plane P r projecting parallel to the Z-axis into the XF-plane. 
We say that P r is maximal if it intersects to and if it projects to the maximal line. 

The maximal line K intersects the projection Pui in projected line segments which 
lifts to line segments and isolated vertices inwfl P r . 

These definitions can now be used to stretch a polygon in a consistent way without 
creating self-intersections. 

3.2.1. Stretching a polygon in its maximal line: Suppose L is a line given by x + y = I 
(with I e N) which intersects the projected polygon Put. The maximal line K has 
equation x + y = k, and necessarily, k > I. 

We say BCC basis vectors Cj are transverse to the plane P r if they are not parallel 
to P r . Let Ci be transverse to P r — then e {c 1; c 2 , c 5 , c 6 }; the other lattice vectors 
lie in P r . Each line segment in this intersection can be translated one step in the Cj 
direction by BCC elementary moves using the construction in Figure 6. We restrict 
such translations to be in the Cj direction (with i = 1 or i = 2) in what follows. This 
will translate all line segments of non-zero length in the plane P r in the direction. 

This leaves the case of isolated vertices in P r f]uj. The edges incident to and from 
these vertices will be of the form C\C§ or C2C5 (since K is maximal). In this case one 
may perform the move CiC6 —> C1C1C6C5 or C2C5 — > C2C1C5C6 to translate the isolated 
vertex in the Ci direction. A similar construction will translate the segment in the c 2 
direction instead. 

Observe that no self-intersections can occur since all parts of the polygon in F r nw 
are translated in parallel in the Ci direction (with i = 1 or i = 2). 

Once all line segments in the maximal plane P r are translated in the ci (or 
c 2 ) direction, the polygon lj is said to have been stretched in the plane P r in the 
direction ci. 

On completion of this construction, the maximal line K is translated to the line 
L by translating the plane P r in the — Cj direction. This plane is denoted by Q r and is 
parallel to P 2 . Then Q r projects to L with formally x+y = £ = k — lin the XF-plane. 
L lifts to Q r , and Q r n uj is a collection of line segments and isolated vertices of lo. 

Observe that the departing edges from Q r to the maximal side arriving in P r are 
all in the c» direction, since the polygon was stretched in that direction in the previous 
step. 

In addition, the P r contains no line segments, since these were translated into 
the Ci direction. Furthermore, edges incident arriving in P r from the opposite of the 
maximal side are parallel or anti-parallel to Cj with j = 1 or j = 2. 
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3.2.2. Stretching a BCC polygon We proceed by recursively executing the stretching 
of u> in the Cj direction in the plane Q r . The intersection Q r n ui is a collection of line 
segments and isolated vertices in w. 

If Q r = P r , where P r projects to the maximal line K, then the situation is as 
described in the last section: We must consider two different cases in translating parts 
of the polygon in P r in the Cj direction. The first case involves line segments in P r r\uj 7 
the second case involves isolated vertices in P r n ui. Such a isolated vertices must have 
incident edges CiC & or c 2 c 5 (so as not to collide with P r ). These two cases were already 
dealt with above. 

In the event that Q r ^ P r , suppose that the stretching was recursively done 
starting in P r and moving the plane in the — Cj direction so that the last stretching 
was done in the Cj direction in the plane Q' r = Q r + Cj. Without loss of generality, 
one may suppose that the stretching is done in the i = 1 direction. Then all edges 
between Q r and Q' r are parallel or anti-parallel to ci. 

There are three different cases to consider. The first case involves line segments 
in Q r C\uj, and the second case involves isolated vertices in Q r n cj with incident edges 
C1C6 or C2C5. These two cases are done by translating the vertices and edges in the Ci 
direction similarly to those line segments and isolated vertices in the plane P r above: 
Since line segments projects to lines in the square lattice and BCC moves to BFACF 
moves in the square lattice, these line segments can be translated by applying BCC 
moves in the ci direction. Translating isolated vertices in the second case is similarly 
done. 

The third and final case involves isolated vertices v of the polygon in Q r with 
edges on either side of Q r . Suppose that these vertices are to be moved in the Ci 
direction. Then the arrangements must be civci, or c 2 vci; the middle vertex v in 
these cases lies in Q ri and the second edge moves from Q r to a plane Q' r = Q r + ci. 

In the cases civci the edges are left unchanged, since the departing edges to Q' r 
are already in the ci direction. The case c 2 vci is updated to Ciwc 2 , with w in the 
plane Q' r . Observe that w is always not occupied in Q' r before the move, because all 
arriving edges from Q r to Q' r are in the Ci direction before the moves are done. 

In each case, these constructions give a polygon with departing edges from Q r 
to Q' r in the Cj direction. Observe that the relative orientation of edges are again 
maintained, and that the moves are possible without any self-intersections in u. 

The implementation of the stretching is recursive. Start in the plane Q r = P r 
which projects to the maximal line, and stretch the polygon a number of times in the 
Cj direction transverse to Q r . Then define Q' r — Q r and Q r — > Q r — c, recursively. 
This will stretch the polygon any desired length in each plane Q r intersecting parallel 
to P r without creating an intersection in the polygon. 

The effect of the construction is to cut lj along a plane Q r and to move the 
two parts of lj on cither side of the polygon any number of steps in the c, direction 
transverse to Q r apart while inserting edges parallel or anti-parallel to Cj to reconnect 
it into a single polygon. 

Observe that the construction does not change the knot type of co, and that it is 
the realisation of an ambient isotopy on the complementary space of the polygon. We 
say that we have stretched the polygon u in the Cj direction along the plane Q r . This 
completes the construction. 

Observe that a similar construction will enable one to stretch ui in directions 
transverse to lattice planes in the BCC transverse to any of the basis vectors Cj of 
the BCC. This follows by exchanging the lattice basis vectors, such that a rotation of 
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Figure 7. Stretching a projection of a BCC polygon. In this illustration, a 
polygon was stretched recursively in the projected directions until each projected 
edge was doubled. 

the polygon preserving its chirality, is done in the above analysis. These observations 
complete the proof of the following theorem: 

Theorem 2 Let lo be an (unrooted) polygon in the BCC lattice with projection Puo in 
the XY- YZ- or XZ -plane. Let L be a line in the AB-plane with formula A± B = £ 
where (A, B) is one of(X,Y), (Y,Z) or (X,Z). 

Suppose L intersects the projection Put and that £ is an integer. Suppose also that 
L lifts to the plane Q r parallel to the normal of the projection plane. Then lo can be 
stretched along the plane Q r in a direction Cj transverse to Q r by performing the BCC 
elementary moves in Figures 4 and 5 on to. (} 

In Figure 7 we illustrate that a polygon can be recursively stretched until each 
projected edge has been replaced by two edges in the same projected direction. This 
is called the subdivision of the projected image of the polygon. 

3.3. Contact free polygons 

In this section we show that a given BCC polygon can be swept onto a sublattice 
L with basis vectors {ci, c 3 , C4, c 5 , C7, c 8 }. This sublattice is isotopic to the simple 
cubic lattice, and once the polygon is contained in it, then a subset of the BCC moves 
reduces the usual BFACF moves. 

In order to sweep a polygon to into L, it is necessary to avoid self- intersections 
in the polygon when executing the necessary elementary moves. Such possible self- 
intersections are avoided by stretching the polygon, using the recursive subdivision 
explained in the previous section. 

Two vertices v and w, non-adjacent in a polygon lo (ie not connected by an edge 
in lo), form a contact if v — w = Cj for some value of j. 

Such a contact is said to be a contact in the direction Cj. A polygon lo is contact 
free if it has no contacts. Observe that every contact is a lattice edge and projects to 
a line segment in the XY"-plane. We show that by subdividing a polygon using the 
constructions outlined in the previous sections, one can make it contact free. 

Lemma 3 By applying the BCC elementary moves in Figures 4 and 5 to an unrooted 
polygon in the BCC lattice, it can be transformed into a contact free polygon. 

Proof: Use BCC elementary moves to remove all contacts in the polygon ui as follows. 

If co has contacts in the Cj direction, then subdivide the polygon such that the 
stretching is done in the direction. Since all edges in the Cj direction will be doubled 
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up (effectively the scale of the polygons is changed such that every edge is replaced by 
two edges in a line segment of length 2), this removes all contacts in the Cj direction. 
If there are still contacts in the Cj direction (for some other j), then subdivide the 
polygon in the Cj direction. 

Observe that if there no contacts in the Cj direction, then stretching the polygon 
in the cj direction cannot create new contacts in the c, direction. 

This follows because the creation of a contact in the c t direction will require the 
translation of parts of the polygon relative to one another in the Cj direction, which 
does not occur when stretching in the Cj direction. 

Repeat the process of subdivision, until all contacts are removed and a contact 
free polygon is obtained. <0 

3.4- Pushing contact free BCC polygons into a simple cubic sublattice 

Our goal is to show that every polygon in the BCC lattice can be changed into a 
polygon in sublattice L which is isotopic to the simple cubic lattice by the applying 
the BCC elementary moves. 

We defined L to be that sublattice of the BCC with basis vectors 
{ci, C3, C4, C5, C7, eg}. L is a non-orthogonal lattice, and if lu is polygon in L, then 
the subset of elementary moves in the BCC lattice which only involve the edges in the 
basis of L reduces to the usual simple cubic lattice moves. 

By lemma 3 one can show that BCC polygons can be made contact free and so 
we only have to examine contact free polygons in this section. 

If cj is a contact free polygon in the BCC, then it is pushed into L by removing 
from it all edges in the c 2 or c 6 directions. The method of proof is as follows: Replace 
very edge in the c 2 direction by the three edges C1C4C7 and every edge in the c 6 
direction by three edges c 5 c 8 c 3 . 

All that remains is to show that one can arrange matters such that self- 
intersections will not occur. 

Theorem 4 By applying the classes of BCC elementary moves in Figures 4 and 5 to 
unrooted polygons in the BCC lattice, any such polygon can be swept into a polygon 
in the sublattice L. 

Proof: Let w be a contact free polygon and subdivide it twice in each of the ci , c 5 , 
C3 and C7 directions. Then the shortest distance in the ci direction between vertices 
u, v not connected by an edge of the polygon is at least 3 steps. Similarly for the C5, 
c 3 and C7 directions. 

We now show that edges in line segments in the c 2 or c 6 directions can be swept 
into the sublattice L. 

Replace edges in the c 2 or c 6 directions as follows: 

• If c 2 c 2 c 2 . . . c 2 is a sequence of consecutive edges in the c 2 direction, then replace 
them by c 1 C4C7C 1 c 4 C7 . . . C!C 4 C7. 

• Similarly, replace any sequence of consecutive edges CqCqCq . . . by the sequence 

C5C 8 C 3 C5C 8 C 3 . . . C 5 C 8 C 3 . 

These changes can be achieved by positive non-planar BCC elementary moves on 
the BCC polygon. Observe that these substitutions insert new vertices in the Ci, C5 
and C7 and c 3 directions, adjacent to existing vertices in u>, but that no contacts can 
be created since other vertices in these directions are a distance of at least three away. 
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The only intersections that can arise in the above construction do so at the 
endpoints of a sequence c 2 c 2 c 2 . . . c 2 or c 6 c 6 c 6 . . . c 6 . These are avoided as follows: 

• If the edges at the beginning of the sequence c 2 c 2 c 2 . . . c 2 is Csc 2 . . ., then the 
substitution c 2 — > C1C4C7 on the first c 2 will cause an intersection (or a "spike") 

C5C1C4C7 If instead a neutral non-planar elementary move c 5 c 2 — > C4C7 is 

executed here, then the spike is avoided while the edge in the c 2 direction is 
removed. 

• A similar argument holds if the sequence of c 2 's ends as . . . c 2 c 3 , in which case 
one finds c 2 c 3 — > C1C4. 

• Similar arguments can be used to deal with the case of line segments of the form 
c 6 c 6 . . . c 6 . 

• Lastly, if the line-segment has length one and is of the form c 5 c 2 c 3 , then one 
obtains the negative BCC elementary move c 5 c 2 c 3 — » c 5 CiC4C7C 3 — > C4. 

• A similar argument holds for segments of the form Cic 6 c 7 . 

Completion of these elementary moves produces a polygon with no edges in the 
c 2 or c 6 direction. This is exactly a polygon in the sublattice L. This completes the 
proof. <0 



3.5. Irreducibility classes of the BCC elementary moves 

By Theorem 4 all BCC lattice polygons can be changed into lattice polygons in the 
lattice L which is isotopic to the simple cubic lattice. Since this process is reversible, 
one only has to consider polygons in L and the effects of the BCC elementary moves 
on these polygons in this sublattice. 

Restricting the BCC elementary moves to the sublattice L implies that all 
elementary moves including either the directions c 2 or c@ must be excluded. Observe 
that if both c 2 or c 6 are excluded from the set of BCC elementary moves (for example, 
ci <-» c 2 cic 6 or ci <-» c 2 c 3 c 8 , are excluded) then the remaining moves reduce to the 
standard simple cubic lattice BFACF moves in Figure 1 in L. 

Thus if we restrict ourselves to this subset of possible moves on polygons in L then 
we have the following lemma, which is a direct corollary of Theorem 1 or Theorem 3.11 
in reference [16]. 

Lemma 5 The BCC elementary moves, restricted to the sublattice L, applied to 
unrooted polygons in L, have irreducibility classes which coincides with the knot types 
of the polygons as piecewise linear embeddings in R 3 . <0 

Since every unrooted polygon in the BCC lattice can be (rcversibly) reduced to a 
polygon in L, the following theorem is an immediate corollary of Theorem 4 and the 
last lemma. 

Theorem 6 The irreducibility classes of the BCC elementary moves, applied to 
unrooted polygons in the BCC lattice, coincides with the knot types of the polygons 
as piecewise linear embeddings in M 3 . <0 

This completes the proof. In other words, it follows that the set of BCC 
elementary moves, applied to unrooted polygons in the BCC lattice, is irreducible 
within the knot type of the polygon. 
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Figure 8. The four plaqucttes adjacent to an edge in the FCC lattice. 



4. BFACF-Style Elementary Moves in the FCC Lattice 

In this section we introduce local elementary moves of a BFACF-style algorithm on 
the FCC lattice; in fact there is only a single move that replaces two edges by a single 
edge and vice-versa. As was done for the BCC moves described above, we show that 
this new move is sufficient to realise a piecewise linear orientation preserving isotopy 
between two unrooted FCC polygons of the same knot type. This demonstrates that 
the irreducibility classes of the FCC elementary move coincide with the knot types of 
unrooted FCC polygons as determined by their embeddings in three space. 

Local elementary moves in the FCC lattice have been explored previously in the 
literature, see for example references [7], however, the ergodicity properties of these 
elementary moves in the FCC have not been studied for knotted polygons. 

4-1- The FCC elementary move 

Vertices of the face centered cubic (FCC) lattice are points in K 3 with positions given 
by the linear combinations pdi + qd- 2 + rd 3 + sd^ + td 5 + udg where p, q, r, s,t,u G Z, 
and where the vectors dj is given 



We have labelled these vectors so that d^+g = — d^. We define adjacent, lattice edges, 
end-vertices, lattice polygon and line segment on the FCC lattice in the same way that 
we did on the BCC lattice. Note that a polygon of n edges on the FCC has geometric 
length 

Observe that the sublattice generated by any set of three different non-coplanar 
vectors in the generating set of the FCC is isotopic to the simple cubic lattice. For 
example, the set {di,d 3 ,d5} generates a sublattice of the FCC which is isotopic to 
the cubic lattice (and may be viewed as a non-orthogonal simple cubic lattice). 



di = 
d 4 
d 7 = 

dio 



(1,1,0), 
(1,0,-1), 

(-1,-1,0) 
: (-1,0,1), 



d 2 : 
d 5 
d 8 : 
dn 



(1,-1,0), 
(0,1,1), 

(-1,1,0), 

= (0,-1,-1) 



da = (1,0,1), 
d 6 = (0,1,-1), 
d 9 = (-1,0,-1), 

di2 = (0,-1,1) 



15 



Figure 9. Translating a line segment in the d; direction in a polygon in the dj 
direction using the elementary move in the FCC lattice. 

Similarly, there are sets of three coplanar vectors which generate a two- 
dimensional triangular lattice; for example, the set {di, d 3 , di 2 } are the basis vectors 
of a triangular lattice in a plane with normal vector di x d 3 (and observe that 
di2 = d 3 di). 
FCC Elementary Move: 

• See Figure 8. Any vector di in the FCC lattice is incident to four triangular 
plaquettes. That is, there are four different pairs of vectors (dj, d k ) so that 
di = dj +d k . The substitution of di by djd k is a positive atmospheric FCC move, 
or a positive FCC elementary move and it increases the length of the polygon 
by one edge. Similarly the reverse move, replacing djd^ by di is a negative 
atmospheric FCC move or a negative FCC elementary move and it decreases the 
length of the polygon by one. 

4-2. Stretching polygons by using the FCC elementary move 

In examining the irreducibility classes of the FCC elementary move, we shall follow 
the same strategy used for the BCC lattice in Section 3. 

We first show that FCC polygons can be stretched, then made contact free, and 
finally swept into a simple cubic sublattice of the FCC. In this simple cubic sublatticc 
of the FCC, the usual simple cubic lattice BFACF moves can be performed on the 
polygon by using combinations of the FCC elementary moves. Theorem 1 can then 
be used to complete the proof that the irreducibility classes of the elementary move 
coincides with the knot types of the polygon. 

The basic construction in the proof is illustrated in Figure 9. By using the 
elementary move, the construction in Figure 9 shows that one can translate an entire 
line segment (in the di direction one step in the dj direction, provided that dj ^ ±dj . 

This construction can be performed regardless of the orientation of the edges 
incident at the ends of the line segments, provided that the set of target vertices in 
the dj direction are not occupied by the polygon. 

Let {d i; dj, dfc} be a triple which form a triangular plaquette in the FCC lattice 
(so dfe = dj — dj. Note that these three vectors generate a triangular sublattice T in 
the FCC; this sublattice is a geometric plane A with normal di x dj. 

Consider the (non-orthogonal) projection of an FCC polygon uj into the plane 
A along a lattice direction d; which is transverse to the plane A (that is, d; is not 
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Figure 10. The vectors {di , dj , dj,} generate a triangular lattice in the plane 
A. The Cartesian coordinate system (X, Y) is set up with X-dircction given by 
di + dj and Y-direction given by d^. Observe that |d n | = \/2, so that k = \/3/2 
and I = 1/V2. 

coplanar with vectors {dj, dj, dfc}). We say that this projection is taken along d; into 
the djdj -plane A. 

In general, the projection of an FCC lattice polygon lj along d; into the plane A 
is a lattice knot projection into a triangular lattice which is the projection of the FCC 
onto A. Define a two dimensional cartesian coordinate system in the plane A with 
Y-direction given by d^ and Y-dircction by dj + dj. This is illustrated in Figure 10. 

If Put is the projection of the polygon cj in the plane A along a transverse lattice 
direction d;, then Put is a finite graph in the triangular lattice in A. 

Hence, there exists a right-most line K parallel to the Y-axis (see Figure 10) 
which intersects the projection Put. We say that the line K is the maximal or right- 
most line cutting the projection Put. The right-most line is itself the projection of a 
plane P r projected with co along the (non-orthogocal) directopm d ; . The intersection 
of this plane with the FCC lattice is a triangular sublattice of the FCC, generated by 
the set of basis vectors {d^, d;}. 

Line segments in co are projected to line segments or to points in K, and these 
projected images lift back up to parts of the polygon in the intersection P r n uj of 
the plane and the polygon. Each line segment in this intersection can be translated 
one step in the dj (or dj) direction (transverse to the plane P r ) by using the basic 
constructions in Figure 9. If all the line segments in the plane which projects to the 
maximal or right-most line K are moved in the dj direction, then we say that the 
polygon is stretched in the dj direction from the plane P r . 

Observe that once a polygon has been stretched in the plane P r , then there are no 
line segments in the polygon contained in P r . This, in particular, implies that P r Out 
is a collection of isolated vertices in w. At each such vertex w either passes through 
the plane to its right in the dj direction, or it turns to stay to the left of P r . Projected 
images are illustrated in Figure 11. 

Next we stretch the polygon in a plane P r recursively starting in its right-most 
line K, and then successively moving left. Line segments in P r n lo can be moved 
in the dj direction using the construction in Figure 9. Since all the edges incident 
with P r on the maximal side of P r are in the dj direction, these constructions can be 
performed without creating any self-intersection. 
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Figure 11. The projection of the polygon close to the maximal or right-most line 
K during the stretching of the polygon in the dj direction. The conformations 
of edges incident with the plane P r arc in general one of the cases above. In two 
cases the polygon passes through the P r . In the third case it touches it in one 
vertex before turning back. In this last case the vectors d n and d; must make 
an angle of 60° with one another. Cases as in (a) are removed as illustrated, 
while case (c) is left unchanged. In case (b) the edge in the d„ direction is passed 
through the plane P r in the direction dj. Sec Figure 12 for more details. 

This leaves the cases of isolated vertices in P r n to. The situation is as illustrated 
in Figure 11: All the parts of the polygon to the right of P r K been stretched in the 
dj direction, and the next step is to move parts of the polygon which has isolated 
vertices in P r in the dj direction. 

In Figure 11(a) we depict a situation in which the two vectors a„ and a/ must 
make a 60° angle with one another at the vertex in P r . In this case we can remove 
these two edges from the polygon by making a negative atmospheric move replacing 
them with a single edge as shown. Thus we can eliminate this situation and we do 
not need to consider it in the discussion below. 

Next, we address the situation depicted in Figure 11(b). We need to translate 
the edge d n in the dj direction. This is done as illustrated in Figure 12 (see also 
Figure 11). If the vertex marked by o in the left-most figure in Figure 12 is vacant 
(that is, not occupied by the polygon), then we can use two elementary moves to 
change the conformation d n dj to djd„; so the vector d„ is translated in the dj 
direction, as desired. 

On the other hand, if the vertex marked by o is occupied by another part of 
the polygon, then there are only two possibilities. First, it may correspond to the 
conformation depicted in Figure 11(a), in which case the elementary move can be 
used to remove the occupied vertex. 

Otherwise it is the situation shown in the middle of Figure 12. In this case let Q 
be the plane coplanar with {dj,d„} and consider the line S — Q n P r . Move along 
S in P r until an open vertex o is enountered as shown the second and third parts of 
Figure 12. 
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Figure 12. The case in Figure 11(b). If Q is a plane coplanar with {dj, d n } 
then S is the line S = Q D P r - If one moves along S in the vertical sense here, 
then eventually one must encounter an open vertex o. In this case the moves 
in Figure! 11(b) can be systematically performed on each of the isolated vertices 
along S, starting at the top. This will translate all the edges on the left of S 
through the plane P r towards the right in the dj direction. 

Since the polygon is finite we eventually must find a vacant vertex — as shown 
in the rightmost past of Figure 12. The edge (depicted as dj) can then be translated 
in the dj direction. This (effectively) moves the vacant vertex back along P r and so 
one recursively apply this construction until we finally move the desired edge, d„, in 
the dj direction, using the construction in Figure 11(b). 

The stretching of the polygon in a new plane, Q ri in the transverse (to Q r ) 
direction dj, proceeds by first finding the right most line in the projection, and lifting 
it to the plane P r (parallel to Q r ). The polygon is then stretched one step in the dj 
direction in P r . Then P r is moved closer to Q r one step in the — dj direction and the 
polygon is stretched recursively until it is finally P r is coincident with Q r in which 
case the polygon is stretched one step in the dj direction on one side of Q r . 

The effect of this construction is to cut the polygon along the plane Q r into a left 
part and right part. The right part is then translated in the dj direction and the two 
parts are reconnected by inserting inserting edges in the dj direction. 

In topological terms this construction is an (orientation preserving) ambient 
isotopy of three space, and it does not change the knot type of the polygon. We 
say that the polygon u> was stretched in the dj direction transverse to the plane Q r . 

This construction proves the following lemma: 

Lemma 7 Let lj be any FCC lattice polygon and let Q r be a lattice plane intersecting 
iv. Then uo can be stretched, using the FCC elementary move, in Q r in the dj direction 
if dj is transverse to Q r . 

Proof: Since Q r is a lattice plane, determine a plane A transverse to it which intersects 
the FCC in a triangular lattice. Orient he projected polygon as above, and then use 
the elementary moves to stretch the polygon in Q r in the desired direction as described 
by the construction above. <0 

4-3. Contact free FCC lattice polygons 

As was the case for the BCC lattice, care must be taken to ensure that different pieces 
of the polygon do not come too close together while we try to sweep it onto a simple 
cubic sublattice. Thus we define contacts on the FCC lattice in the same way we did 
on the BCC lattice; two non- adjacent vertices u and v in a polygon u> form a contact 
if u v = dj for some j. That is, these vertices are a distance \f2 apart in the FCC. 
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A polygon uj is contact free if it has no contacts. We next show that we can use 
the FCC elementary move to transform any polygon into a contact free polygon. 

Lemma 7 shows that any given polyon can be stretched in the dj direction 
transverse to a plane Q r (which intersects the FCC lattice in a triangular sublattice) . 

If a polygon is stretched in every transverse direction in every lattice plane Q 
which intersects it in a triangular lattice, then the effect is to double the length of 
every edge of the polygon. That is, each edge gets replaced by two edges. This is a 
subdivision of the polygon uj. One may repeat this stretching several times to stretch 
each edge any number of times. 

If polygon has a contact in the dj direction, then it is removed if the polygon 
is subdivided once since all distances between vertices in the polygon is increased by 
at least a factor of 2 in a subdivision of the polygon. In addition, a contact in the 
dj direction cannot be created if the polygon is stretched in the dj direction, since 
the relative orientations of vertices and edges in the dj direction is maintained if the 
stretching is in the dj direction. 

Thus, one may remove all contacts from a polygon by stretching or by subdivision. 
This leaves a contact free polygon, and we proved the following lemma. 

Lemma 8 By subdividing an FCC polygon, it is possible to transform it into a contact 
free polygon, using the FCC elementary move. 

Proof: Consider a (non-orthogonal) projection Puj of an FCC polygon uj into a lattice 
plane A along a lattice direction transverse to A, where A intersects the FCC in a 
triangular lattice L t with basis vectors (say) dj, dj and d^. 

The projection Puj cuts the plane A into one infinite face and a set of finite faces 
(or areas). Each finite face has an area measured in units of the elementary triangle 
in the triangular lattice. 

Since dj is a basis vector of the sublattice L t , a subdivision of to in the dj direction 
increases the area of each finite face in the projection in the plane A by at least one 
unit elementary triangle. This is similarly true when subdividing in the dj and d^ 
directions. 

Contacts in the FCC polygon which are oriented in the dj, dj and d fc directions 
will project in the plane A as edges in L t . A subdivision in the direction of a given 
contact will remove it from the projection, and thus also from the polygon itself. 

If subdivisions in the dj direction removed all contacts in this direction, then 
subsequent subdivision in other directions will not create contacts in the dj direction 
(since this would require the translation of parts of the polygon in the dj direction, 
and this cannot occur). 

Hence, by using subdivisions of uj in each of the six independent directions of the 
FCC lattice, all contacts are removed from the polygon, and it becomes contact free. 
This completes the proof. <0 

4-4- Pushing contact free FCC polygons into a simple cubic sublattice 

Let L be that sublattice of the FCC with basis vectors {di, d3, ds, d7, dg, d xl j. Then 
L is ambient isotopic to the simple cubic lattice. 

If uj is a polygon in L, then the (simple cubic lattice) BFACF move (see 
Figure 1) can be performed on uj in L by composing two FCC elementary moves 
in order to execute each single BFACF move. For example, the positive BFACF move 
di — > d 3 didg can be performed in the two step sequence di — > d 3 d 6 —¥ d 3 did g . 
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In other words, the irreducibility properties of FCC lattice polygons in the 
sublattice L is determined by the irreducibility properties of the simple cubic lattice 
BFACF moves, and by Theorem 1. We shall use these facts to complete the proof 
that the irreducibility classes of the FCC elementary move on unrooted FCC lattice 
polygons coincide with the knot types of the polygons as piecewise linear embeddings 
in three space. 

The strategy is as follows: If lj is an FCC lattice polygon, then it will be moved 
to a polygon in the sublattice L, from where it will be put in standard position A 
given its knot-type K. Since the composition of all the FCC elementary moves is a 
continuous piecewise linear orientation preserving isotopy of three space, we know that 
the polygon A in standard position has the same knot type K, and moreover, that 
every polygon can be moved to A using a finite number of FCC elementary moves. 

It only remains to prove that every FCC lattice polygon can be moved to a 
polygon in L. 

Lemma 9 By using the FCC elementary move, any FCC lattice polygon to can be 
swept into a sublattice L of the FCC, where L is ambient isotopic to the simple cubic 
lattice. 

Proof: Since every polygon can be made contact free, it is enough to assume that to is 
contact free. We next prove that every such contact free polygon can be moved into 
the lattice L. 

If lo is embedded in L, then we are done. 

Otherwise, to has some edges in directions from the set {d 2 , d 4 , d 6 , d 8i di ,d 12 }. 
Since w is contact free, the FCC elementary move can be performed on each of each 
of the edges in w without creating self-intersections. 

Proceed along the polygon and remove edges in the dio d 2 , d4, d 6 , d 8 , dio and 
di2 directions systematically by replacing them as follows: 

d 2 -> d 3 dn, d 4 -> didn, d 6 -> dgd 5 , 
d 8 -> d 9 d 5 , dio -> d7d 5 , d i2 -> d 3 d 7 . 

If during this process a contact is created, then the polygon is again swept in the 
{di,d3,d5} directions. Each contact will have a non-zero projection along at least 
one of these directions, and will then be removed without inserting any new edges 
in any directions apart from those in the simple cubic sublattice L. Finally, these 
elementary moves will move ui into the sublattive L. Q 
A direct corollary of this lemma is 

Theorem 10 The irreducibility classes of the FCC elementary move, applied to 
unrooted polygons in the FCC, coincides with the knot types of the polygons as piecewise 
linear embeddings in M 3 . 

Proof: Observe that since L is a sublattice of the FCC which is isotopic to the 
simple cubic lattice, the irreducibility classes of the FCC elementary move applied 
to unrooted polygons in C coincides with the knot types of the polygons as piecewise 
linear embeddinds in M 3 . This follows directly from the observation above that all 
BFACF moves on a polygon in L can be induced by using the FCC elementary move, 
and from Theorem 1 or Theorem 3.11 in reference [16]. 

By Lemma 9 any unrooted polygon in the FCC can be (reversibly) swept into the 
sublattice L using the FCC elementary move. 
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Standard BFACF elementary moves in the sublattic L can be performed by 
combining FCC elementary moves. For exampl, the positive move di — > dsdidg in L is 
a BFACF move, and it is performed in two steps in the FCC by di — > d 3 di 2 — > d 3 did 9 , 
since di = d 3 + di 2 and d i2 = di + d 9 . The other cases, and neutral and negative 
SC BFACF moves can similarly be checked. 

Hence, since the standard SC BFACF moves in Figure 1 can be performed 
on the polygon in L, using composite FCC elementary moves, the theorem follows 
immediately. <0 



5. Conclusions and implementation 

In this paper we have shown that the BFACF elementary moves on simple cubic 
lattice polygons can be generalised to similar moves on the BCC and FCC lattices. 
For the BCC lattice, we propose a set of four elementary moves (two are neutral, and 
two are positive/negative), while on the FCC lattice we propose a single (reversible) 
elementary move (of type positive/negative). 

We proved that, on the BCC and FCC lattices, these moves are suitable to give 
irreducibility classes of unrooted polygons that coincide with the knot types of the 
polygons. This is, we have generalised Theorem 3.11 in reference [16] to both the 
FCC and BCC lattices. For the BCC lattices, our method of proof relied on all four 
proposed elementary moves and so we know that these moves are sufficient. It is 
not clear that any of these are necessary, or that a smaller set of similarly defined 
elementary moves may be sufficient. 

Using the proposed moves we have implemented the GAS algorithm [14] to sample 
knotted polygons on both the FCC and BCC lattices. We give some results in Table 1 
and further numerical results will be published elsewhere. 



5.1. The GAS implementation of the BCC and FCC elementary moves 

In the GAS implementation [14] of the algorithm, lattice polygons of (fixed) knot type 
K are sampled along a Markov Chain by executing the elementary (or atmospheric) 
moves on the polygons. Let <p — (4>n) — (4>o, <fii, <p2> ■ ■ ■) De a realisation of a sequence 
started in the initial state 0o, which is a polygon of length £(<fio) and of knot type K. 

A state 4> n +i from <fr n is generated as follows: Let A + (<j> n ) be the set of all 
positive atmospheric or elementary moves on the polygon (f> n , and similarly, A (4> n ) 
and A_ (</>„) the set of all neutral and negative atmospheric moves on <p n . These are 
the positive, neutral and negative atmospheres of <p n - 

Define the sizes of the of the positive, neutral and negative atmospheres of (f> n by 
a + (4> n ), ao(4> n ) and a_(0„). If £((f> n ) = N, then each of these statistics has a mean 
value denoted by (a+)jv, (o-o)n and (a_)jv- 

Next, define a set of parameters of the GAS algorithm by 

(a_) N 

With these definitions, one may now determine <f> n +i- 

The state 4> n +i is obtained from <p n by selecting with probability P + a positive 
atmospheric move, with probability Pq a neutral atmospheric move, and with 
probability P_ a negative atmospheric move, and executing it on <j> n to obtain (f> n +i- 
In each case, once a determination is made to execute a move, the move is selected 
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uniformly from amongst the available moves in each of the sets A + (<f> n ), A ((f>) and 
A-(<t> n ). 

The probabilities P + , P and P_ are given by 

p = P N a + ((j} n ) 

+ a-(<j>„) +a (<j> n ) + (3 N a + {(f> n y 
p _ ao{(f>n) . 

a-(<p n ) +a (<p n ) + l3 N a + (4> n y 

p = ai{4>n) 

a-(<j>n) + aoOn) + /3jva+(0„)' 
where N — £(4> n ) is the length of the polygon </>„, and (3n is estimated by updating it 
recursively in the simulation. 

This implementation give a sequence (p = {4>)m of M + 1 states sampled from 
polygons of knot type K. The analysis of the data sequence depends on a weight W((f>) 
of the sequence, which is defined as follows: Define a(n,n + 1) = — 1 if <j> n — »■ 
via a positive atmospheric move, and a(n, n + 1) = 1 if <j> n — > <p n +i via a negative 
atmospheric move. Then the weight of the sequence <p is given by 



a_((t) ) + a (<f> ) + (3 e(M a + (<j) ) 1 



n,n+l) /o\ 

(0) " [«-(^)+ao(0L) + ^)a + (^)J 11 V.) ' 

where is the length (number of edges) in polygon <^ n , and where L = \4>\ is the 

length (number of states) in the sequence <j). 

The expected value of the weight of sequences ending in states of length N is 
(W(4>)}n, and the basic GAS theorem states that 
E\ t \=n(W(t)) n = PN (K) 

E\ p1 =m(W( P ))m Pm(KY [ ' 

where r are all possible sequences of states ending in a state of length |r| = N, and p 
are all the possible sequences of states ending in a state of length \p\ — M. 

In this implementation, the FCC and BCC elementary moves described in this 
paper will give an approximate enumeration algorithm, and estimates of the number 
of polygons of knot type K and length n, p n (K), can be obtained. This approach 
produced the data in table 1. 



5.2. The BFACF-style implementation of the BCC and FCC elementary moves 

The GAS algorithm is by no means a standard simulation method (though the authors 
hope that it might become more popular) , and so we finish the paper by discussing how 
the proposed moves may be implemented in a BFACF-style algorithm, as implemented 
in reference [2, 5]; further implementations can be found in references [6, 18]. 

In general the BFACF algorithm is applied to either a self-avoiding walk or self- 
avoiding polygons and so we will discuss the implementation of the BCC and FCC 
moves to sample either walks or polygons. Since these moves do not translate the 
endpoints of a self-avoiding walk, assume that the walk is rooted at its endpoints. 
The simple cubic lattice moves are not irreducible on the state space S$ x of walks 
with fixed endpoints (the origin) and a lattice site x [8, 16] unless those endpoints 
are a distance of at least 2 apart (for details, see reference [8]). We believe that a 
similar result will hold for the BCC and FCC moves we have proposed, but we have 
not yet investigated this. 
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Let uj n be the current walk or polygon (or "state") of length \u n \ edges. Choose 
an edge s from w„ with uniform probability. Enumerate the possible elementary moves 
on s; in the BCC lattice there are 12 possible moves while in the FCC there are 4 
possible moves. At most one these elementary moves is a negative atmospheric move 
(will shorten the polygon or walk), while the remaining moves are positive or neutral. 

Choose one of the possible moves so that a particular positive move is done with 
probability P + , a particular negative move is carried out with probability P_ and a 
particular neutral elementary move is performed with probability P . This produces 
a state uj' n . If this state is self-avoiding then we accept it and set to n +i = co' n and 
otherwise we reject it and set uj u +i = u n . It remains to specify the probabilities P + , 
P- and Po, and hence a probability is attached to each of these moves. This must be 
done differently on each lattice. 

Consider a polygon or walk in the BCC lattice and a given edge s. The possible 
moves we can perform on s depend on the conformation of the of s and its predecessor 
and successor edges. In particular, the moves on a given edge are in one the following 
classes: 

• All the possible moves on the edge s are positive elementary moves increasing the 
length of the walk or polygon; this gives 12 positive elementary moves on s. 

• The possible elementary moves include 11 positive and 1 neutral move; denote 
this as (11, 1,0). 

• The possible elementary moves include 10 positive and 2 neutral moves; denoted 
by (10,2,0). 

• The possible elementary moves include 9 positive, 2 neutral and 1 negative moves; 
denoted by (9,2, 1). 

• Similarly, the following combinations are possible (9,3,0), (8,4,0), (7,4,1) and 
(6,6,0). 

In addition, one has the Boltzmann type relation P + = e~ 2 "P_ relating the 
probability of making a positive elementary move, to that of making a negative 
elementary move. The cases above sets up a set of bounds on the different probabilities: 

11P++P) <i; 

9P++3P <i; 

8P+ + 4P <1; 

6P++6P <1. 



12P+ < 1 

10P+ + 2P < 1 

9P+ + 2P + P_ < 1 

7P+ + 4P + P_ < 1 



We need to maximise the value of P+ while satisfying these bounds . Since f3 > and 
P + = e~ 2/3 P_, it follows that P + < P_. If we additionally assume that P + < P < P_, 
then the first six inequalities are redundant and one can simply study the last two 
(that is, the bottom row). Maximising then gives 

e -2/3 l + e -2/3 1 

P+ = 3(l + 3e- 2 0)' P ° = 6(1 + 3e- 2 ^)' P ~ = 3(1 + 3e-^)' (5) 

Note that P = \{P+ + P_). 

The situation is simpler on the FCC lattice due to the simplicity of the elementary 
move. When an edge s is selected from the polygon there are only two possibilities: 

• the four possible moves on the edge s are positive elementary moves 

• three of the possible elementary moves are positive and one is a negative 
elementary move. 
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This gives the following set of inequalities for the probabilities P + , and P_ : 

4P+ < 1; 3P+ + Pi<l. 

Since there are no neutral moves in the FCC algorithm we do not require Pq (or we 
can simply set Pq = 0). Again we have a Boltzmann factor relating the probabilities, 
namely P + = e~Pp_ with (3 > 0. Since P + < P_ the first of these equations is 
redundant, and the second gives 

e~ p 1 
P+ = l + 3e-20' P - = l + Ze-P- (6) 

With the above choices of {P+, Pq, P-} on cither the BCC or FCC lattices, the 
algorithm for polygons runs as described below. The algorithm for walks is identical 
after the initialisation step. 

1. Initialise the algorithm by choosing the first polygon, u>o, to be a polygon of given 
length and desired knot type K. Define the parameter (3. The length of this 
initial polygon is not crucial except that it must be sufficiently long to embed a 
knot of type K. 

2. Let co n be the current polygon of length |u> n | edges. With uniform probability, 
l/|ix)„|, choose an edge s in tj n . Enumerate all possible elementary moves which 
can be performed on that edge. 

3. With probability P + perform a (particular) positive atmospheric move, uniformly 
selected from amongst all the possible positive moves. Similarly, with respectively 
probabilities Po or P_, perform a particular neutral or negative move. If the sum 
of these probabilities is less than 1, then the remaining default move is to leave 
the polygon unchanged. 

4. Step 3 induces a proposed polygon oj' n . If the proposed polygon is not self- 
avoiding, or if (by default) no elementary move occurred, then uj' n = uj n . 

5. Put w„+i = co' n to find the next state, and continue at step 2 above. 

In this implementation the transition probability of obtaining a polygon r from 
a given polygon cj, P^ T , is given by 

' [P + /|w|], if the elementary move is positive; 

P^^t = \ [-P-/M]: if the elementary move is negative; (7) 

^ [Po/M], if the elementary move is neutral, 

since the edge S is selected with a priori probability l/|w| before the elementary move 
is performed. These transition probabilities satisfy the condition of detailed balance 
between states to and r given by 

M P^ T = \t\P t ^, (8) 

with the result that, in the BFACF-style implementation, the algorithm will sample 
asymptotically from the distribution 

where p n {K) is the number of polygons of length n edges and knot type K, since the 
algorithm is ergodic on the class of unrooted polygons with fixed knot type K. This 
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is not the canonical Boltzmann distribution over the state space of unrooted polygons 
with knot type K given by 

Note that denominator is finite when j3 < — log fix where [Ik is the growth constant 
of polygons of given knot type defined by 

limsup [p n (K)] 1/n = fi K . (11) 

This limit is known to exist for the unknot, but its existence for other knot types is a 
long standing open problem. 



5.3. Metropolis- style implementation of BCC and FCC elementary moves 

A simpler implementation of the BCC and FCC elementary moves uses the Metropolis 
algorithm. In this case steps 2 and 3 above are modified as follows: Choose an edge 
s with uniform probability l/|w„| in the current polygon, and select (uniformly) one 
of the 12 possible elementary plaquettes on the edge in the BCC lattice, or one of 
4 possible elementary plaquettes on it in the FCC lattice. Then choose one of the 
elementary plaquettes uniformly and attempt the elementary move it defines. If the 
resulting object is self-avoiding polygon oj' n of length \u' n \ edges, then accept it with 
probability min{l, e^'"™'"'"™')} as the next state tJ n +i- 

Otherwise, put u> n +i = w„. In this implementation one can show that the 
condition of detailed balance is given by equation by equation (8), and if /? < — log /x^, 
the algorithm samples asymptotically from the distribution Pp(s) in equation (9). 
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